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ABSTRACT 

For a tree with the given sequence of vertex degrees the spectral radius of its 
terminal distance matrix is shown to be bounded from below by the the average 
row sum of the terminal distance matrix of the, so called, BFS-tree (also known as 
a greedy tree). This lower bound is typically not tight; nevertheless, since spectral 
radius of the terminal distance matrix of BFS-tree is a natural upper bound, the 
numeric simulation shows that relative gap between the upper and the lower bound 
does not exceed 3% (we also make a step towards justifying this fact analytically.) 
Therefore, the conjecture that BFS-tree has the minimum terminal distance spectral 
radius among all trees with the given degree sequence is valid with accuracy at least 
97%. The same technique can be applied to the distance spectral radius of trees, 
which is a more popular topological index. 
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1 Introduction 


We study simple connected graphs. Let V(G) be the vertex set and -E(G) be the edge 
set of an undirected graph G. For any pair of vertices u,v E V (G) let dciu, v) denote 
the distance (the length of the shortest path) between u and v in G. 

The matrix D{G) := {dG{u,v))u,v£V{G) is known as a distance matrix of a graph 


G. Along with the adjacency matrix and the Laplacian matrix 16 , the D{G) and 
related matrices are the most popular objects to study in algebraic graph theory. 

Denote by dciy) the degree of a vertex v E V (G) in the graph G, i.e., the number 
of vertices being adjacent to v. The vertex v E V (G) is said to be pendent if dciv) = 1. 
All other vertices of the graph G are referred to as internal. By W{G) we denote the 
set of all pendent vertices of a graph G, and let M (G) = V (G) \ W (G) be the set of 
its internal vertices. 

The matrix RD{G) = {dG{u,v))u,v&w{G) is typically referred to as the terminal 
distance matrix or the reduced distance matrix of graph G. A tree is a connected 
graph with N vertices and N — 1 edges. Terminal distance matrices of trees are of 
special interest, since a tree can be reconstructed by its terminal distance matrix 
(see 


22 


Concepts based on the distance matrix are intensively employed in the mathe¬ 
matical chemistry. In particular, one of the oldest topological molecular indices, the 
Wiener index, is dehned as one half of the sum of all elements of the distance matrix 
of a graph: 

VI{G)='t Y, daGv), 

u,v&V{G) 

and represents just an example from the large family of distance-based topological 
indices [^[^. The terminal Wiener index of graph G is dehned by analogy as the one 
half of a sum of elements of RD{G): 


TWI{G) = t Y da{u,v). 

u,v&W{G) 

Spectrum-based indices, which are calculated using eigenvalues and eigenvectors of 
various graph matrices, form a yet another family of topological indices , the most 
famous being the Estrada index [^. Balaban et ah suggested the distance spectral 
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radius (the largest eigenvalue of the distance matrix) as a molecular descriptor giving 
rise to the extensive QSPR[^ research and to the studies of mathematical properties 
of the distance spectral radius {DSR). 

In particular, in upper and lower bounds were suggested (and later improved 
by Zhou [^) for DSR of a tree in terms of the tree order and the sum of squares of 
distances between unordered pairs of vertices. Das obtained bounds for DSR for 
bipartite graphs in terms of partition orders and characterized extremal graphs. 


Indulal 11 has shown that 


DSR{G) > -WIiG) 


n 


( 1 ) 


for any connected graph G with n vertices, with equality if and only if G is distance 
regular, i.e., when all row sums Di, ...,D„ of D{G) (also known as distance degrees) 
are equal to each other. He also proved the inequality DSR{G) > \J with 
equality if and only if G is distance regular, and even stronger lower bounds in terms 
of row sums, the second degree distance seguence, and the Wiener index (see 0 


for details). Alternative bounds of this sort were reported later by He et al 12 


In (2^ a simple lower bound has been suggested in terms of graph order and of two 


maximal vertex degrees. Upper and lower bounds of DSR were also obtained in 27 


for bipartite graphs in terms of graph order, diameter, and extremal degrees in two 


partitions. Recent results are summarized in 18 


Definition 1. 15 ,24 For a real seguence x = (xi, Xn), n G N, denote with 




= (xp], ...,X[„]) the seguence, where all components ofx. are arranged in ascending 


order. 


For a graph G dehne its degree seguence as d(G') := (dG('i^))]^gv"(G)- 

Definition 2. A seguence d = {di,...,dn) is called graphic if such a graph G exists 
that d(G) = d (see j^). If G is a tree, we will say that d generates a tree. 

A non-decreasing natural sequence d = (di, ...,dn) is known to generate a tree if 
and only if di -|- ... + dn = 2{n — 1). Let T(d) be the set of all trees with the degree 
sequence d. 

^QSPR = quantitative structure-property relations. 
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An upper bound for DSR{G) has been suggested in 27 in terms of minimum 
degree, second minimum degree and the diameter of graph G. Later in [2,14 


a series 


of improved upper bounds were formulated for DSR{G) in terms of degree sequence 
and distance matrix row sums of G. 


In 23 Zhang introduced the notion of BFS-tree (breadth-hrst-search tree) and 
proved that BFS-tree has the largest spectral radius of the Laplacian matrix among 
all trees with the given degree sequence. The BFS-tree generated by a degree se¬ 
quence d = {di,...,dn) is denoted as and is built in a “top-down” manner 

by sequentially hlling “levels”. We start with a single vertex at Level 1 and connect 
it to dn vertices at Level 2. Then we add dn-i -l- ... -I- d^-d^ — dn vertices to Level 3 
and connect them sequentially “left to right” to dn vertices from Level 2 so that the 
latter have degrees dn-i, ■■■,dn-d„ respectively. In the same manner we add vertices 
to Level 4 and connect them sequentially “left to right” to the vertices from Level 
3 picking the largest unused degree from the sequence d, and so forth until all n 
vertices are added to the tree (see the example in Fig. [^and 23,24 for the detailed 
algorithm). Later BFS-tree (also known as a greedy tree) was shown to minimize the 


Wiener index 21,24 , the terminal Wiener index 19 , and the maximal number of 


subtrees 25 among all trees with the given degree sequence. 



Fig. 1. An example of BFS-tree for the degree sequence (1,..., 1, 2, 3,4,4,5). 

10 times 

A special case of the BFS-tree for the degree sequence (1,..., l,r, A,..., A), where 
1 < r < A, is referred to as the Volkmann tree. It is shown [^ to be extremal with 
respect to the Wiener index and to many other distance-based topological indices 
over all trees of order n and maximum degree A. 

Stevanovic and Ilic conjectured that the Volkmann tree has the minimum 
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DSR among all trees of order n and maximum degree A. Taking into account the 
above considerations, this conjecture can be generalized as follows: 


Conjecture 1. The BFS-tree has the minimum distance spectral radius among all 
trees with the given degree seguence. 

The terminal distance spectral radius (TDSR) is dehned by analogy to DSR as 
the spectral radius of the terminal distance matrix RD{G) of graph G. This index 
is less studied in the literature, yet its behavior is alike that of DSR (at least, for 
trees). Therefore, a conjecture similar to Conjecture can be put for TDSR: 

Conjecture 2. The BFS-tree has the minimum terminal distance spectral radius 
among all trees with the given degree seguence. 


In the present paper we do not prove these conjectures in full. Nevertheless, below 
we show that DSR of the BFS-tree is at least very close to the minimum DSR of a 
tree with the given degree sequence, and TDSR of the BFS-tree is very close to the 
minimum TDSR. 

More formally, we show in Section that DSR of a tree with degree sequence d 
having n pendent vertices is bounded from below by TLB{d), the average row sum of 
the distance matrix of the BFS-tree (or, equivalently, by the value ‘^WI{BFS{d))). 
On the other hand, the minimal value of DSR does not exceeds TUB{d), which is 
dehned as the distance spectral radius of the BFS-tree. That is, if a tree T* minimizes 
DSR over T(d), we show that 

LBid) := -WliBFSid)) < DSRiT*) < UBid) := DSRiBFSid)). (2) 
n 


Similar inequalities trivially hold for the tree T**, which has the minimum terminal 
distance spectral radius: 


TLB{d) := -TWI{BFS{d)) < TDSR{T**) < TUB{d) := TDSR{BFS{d)). (3) 


n 


After some preliminary dehnitions in Section we provide the alternative justi- 
hcation of inequalities ([^ and (|^ in terms of the Wiener index for vertex-weighted 
graphs, the distance-based index hrst introduced in 13 and later studied in [9,10. 
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In Section we show that TLB{d) is very close to Tt/i?(d), and, therefore, 
TDSR{BFS{d)) approximates the extremal value TDSR(T**) well. 

In the conclusion we outline open issues and perspectives. 


2 Main inequalities 


For a symmetric real n x n matrix A denote with Ai(y4),..., A„(A) its eigenvalues 
enumerated in descending order. The biggest eigenvalue of matrix A is called its 
spectral radius. The corresponding eigenvector is called Perron vector and is denoted 
as p{A). Let Sn be the unit sphere in M„, and let be its intersection with the 
non-negative orthant. 

Since, by inequality Q, DSR{G) > ^WI{G) for any connected graph of order 
n, we conclude immediately that for any degree sequence d generating a tree with n 
vertices and any T G T(d) we have 


DSRiT) > min -WIiG). 
Ger(d) n 


From 21,24 we know that the Wiener index is minimized over T(d) by the 
corresponding BFS-tree. On the other hand, since BFS{d) G T(d), the minimum 
value of DSR trivially does not exceed DSR{BFS{d)), and inequalities ([^ follow 
immediately. 

For TDSR an inequality similar to ([^ is easily obtained. Since, by the Raileigh- 
Ritz principle, Ai(74) = max^g 5 ^for any real symmetric nxn matrix A, we 
conclude that for any connected graph G with n pendent vertices 

Xi{RD{G)) > f RD{G) {- 


n 


\\/n 


1 1 

,..., 


> 


> - (1,..., 1) RD{G) (1,..., if = -TWI{G), 
n n 


and inequalities ([^ are obtained similar to ([^. 


3 Huffman Trees and Generalized Wiener Index 


Below we provide an alternative proof of inequalities ([^ and (|^ using recent results 
[9,10 on minimization of the Wiener index for vertex-weighted trees. In this section 
we recall basic dehnitions and theorems following (Tol. 
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A graph G is called vertex-weighted, if each vertex v G V{G) is endowed with a 
non-negative weight Hciy)- 

"as 


The Wiener index for vertex-weighted graphs is dehned in 13 


VWWI{G):=- ^ p,G{u)fiG{v)dG{u,v). 


u,v&V{G) 


Definition 3. Consider a vertex set V with \V\ = n. Let the function /i : 1/ —)■ 
M+ assign a non-negative weight fi{v) to each vertex v & V, while the function d : 
1/ —>■ N assign a natural degree d{v). The tuple {fi,d) is called a generating tuple if 
Ylv&v ^('^) ~ “ ^)- Denote with Jl ;= Yhv&v the total weight of vertex set V. 


Below we use a degree sequence d as a synonym of a degree function d{-) using 
natural correspondence between these concepts. 

Let WT{fi,d) be the set of trees over the vertex set V with vertex weights /i(-) 
and degrees d{-), and let V{fi,d) be the domain of functions of a generating tuple 
(/i, d). Introduce the set W{fi,d) := {w G V{fi,d) : d{w) = 1} of pendent vertices 
and the set M{fi,d) := V{fi,d)\W{p.,d) of internal vertices. 


Definition 4. We will say that in a generating tuple {pi, d) weights are degree- 
monotone, if for any m,m' G M{pi,d) from d{m) < d{m') it follows that pi{m) < 
pi{rn'), and, also, pi{w) > 0 for any w E W{pi,d). 

For a generating tuple (/i, d) the generalized Huffman algorithm builds a tree 
H G WT{pi, d) as follows. 

Setup. Define the vertex set Vi := V{pi,d) and the functions /i^ and d^, which 
endow its vertices with weights p.^{v) := pi{v) and degrees d}-{v) := d{v), v G Vi. We 
start with the empty graph H over the vertex set V{p,, d). 

Steps i = 1, ...,q — 1. Denote with m* the vertex having the least degree among 
the vertices of the least weight in M{p\ d*). Let wi, ...,Wd{mi)-i be the vertices having 
d(mj) — 1 least weights in W{p\d^). Add to H edges wimi, 

Dehne the set I^+i := ..■,Wd{mi)-i} and functions d*’''^(-), endowing 

its elements with weights and degrees as follows: 


/i*+i(y) ;= ^^{v) for V 7^ mi, := p\mi) + -f ... p\wd{mi)-i), 
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d^^^{mi) := 1 . 


( 4 ) 


:= d\v) for v ^ rrii, 


Step q. Consider a vertex nig G M^jjfl^d^). By construction, \M[ijfl,d^)\ = 1, 
\W{ijfl,d^)\ = dirriq). Add to H edges connecting all vertices from W{ijfl,d'^) to nig. 
Finally, set Huiy) := /i(n), v G V{H). 

In general, the Huffman tree is not unique, since there can be more than one set 
of vertices having d{mi) — 1 least weights in d*) at Step i. Denote with d) 

the set of all Huffman trees for the generating tuple {^,d). 

Theorem 1. [[^ If weights are degree-monotone in a generating tuple (/i, d), then a 
tree T minimizes the Wiener index over the set WT (/r, d) of trees with given vertex 
weights and degrees, if and only if T ^ 'H(/r, d). 


If /i(G) := {pI'g{v))v^v(g) is a vector of vertex weights in graph G, VWWI{G) 
reduces to the quadratic form |/r(G)^D(G)/i(G)|^ 

The “classic” Wiener index hFJ is a special case of VWW I for unit weights. It is 


shown in |^[M| that for T G T(d) we have WI{T) > WI{BFS{d)) with equality if 
and only if T ~ Therefore, from Theorem]^ we conclude that 


Note 1. BFS{d) is isomorphic to some Huffman tree for unit vertex weights and the 
same degree sequence d. 


If we dehne the vector of terminal vertices’ weights as w(G) := {tiG{v))„^wiG), the 
terminal Wiener index for vertex-weighted trees is dehned by analogy to the “classic” 
terminal Wiener index as 


TVWWI{G) = ^ IiG{u)f^G{v)dG{u,v) = ^w{GfRD{G)w{G). (5) 

u,v£W{G) 

This index is a special case of VWW I for internal vertices having zero weights, 
and the terminal Wiener index TWI is a special case of TVWWI for terminal vertices 


having unit weights. Similar to the Wiener index, it is shown in 19 that TWI{T) > 
TWI{BFS{d)) for all T G T(d). Hence, from Theorem]^ we see that 


Note 2. BFS{d) is isomorphic to some Huffman tree for unit weights of pendent 
vertices, zero weights of internal vertices, and degree sequence d. 


^To simplify notation we always assume below that vector components and matrix rows go in the 
same order and no confusion arises. 






Note 3. If in a generating tuple (/x, d) internal vertices have zero weights, then weights 
are always degree-monotone in {ia,d) and, by Theorem^ TVWWI{G) is minimized 
over WT(/a,d) with some Huffman tree. 

4 Lower Bounds: Alternative Proofs 

First we estimate from below TDSR{T) := Xi{RD{T)) of a tree with the given degree 
sequence. If a vector fv = (/ii, and a tree degree sequence d = (di, ...,dp^) are 

given, where n < N, denote with {/i, d) a generating tuple obtained by associating 
weights to degrees ”left-to-right” and assigning zero weights to the rest of the vertices. 

Lemma 1 . If a degree seguence d has n elements being egual to unity, then for any 
tree T G T(d) the ineguality holds 

TDSRiT) > 2maxTVWWI{H{ia,d)), (6) 

where H{p,, d) is any Huffman tree from H^ii, d). 

Proof. Obviously, 

TDSR(T) > min Xi(RD(G)). 

Ger(d) 

By the Raileigh-Ritz principle, 

Xi{RD{G)) = max jjl^ RD{G)ia, 

and the maximum is achieved at /i = p{RD{G)). Since any terminal distance matrix 
RD{G) is positive, by Perron-Frobenius theorem, all components of the Perron vector 
p{RD{G)) are positive, and we can limit maximization to Therefore, 

TDSR(T)> min max pf^RD{G)pi. 

G&T{d) ^(is+ 

By the famous minimax ineguality we only decrease the right side by changing 
the order of taking the minimum and the maximum. Therefore, 

TDSRiT) > max min pf^RD{G)fi. 

/iSSn GeT{d) 

From ([^, we have 

TDSRiT) >2 max min TVWWIiG). 

/iSSn GeWT{ii,d) 
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By Note[^ the latter minimum is achieved at some Huffman tree for the generating 
tuple {fJ^yd), and we obtain ([^. □ 

From (§ , the lower bound ( |^ follows immediately, since, by Note[^ BFS-tiee 
is isomorphic to one of Huffman trees for equal pendent vertex weights, and the 
weight vector •••, 1)^ belongs to Nevertheless, below we prove a somewhat 

stronger result, which may be useful in many respects, namely, that the maximum in 
the right-hand side of ([^ is attained when all vertex weights are equal. 

Theorem 2. If a degree sequence d generates a tree with n pendent vertices, then for 
any vector p G and a tree T G >VT(/r, d) the following inequality holds: 

TWI{BFS{d)) > n ■ TVWWI{T). (7) 

We postpone the proof to the Appendix. Note that the reasoning there (namely, 
the Corollary!^ implies that a similar inequality TWI{BFS{d)) > n^TVWWI{T) 
is valid in an even more restricting environment when the positive vector /i of pendent 
vertex weights is taken from the u-dimensional simplex /xi /i„ = 1. 

We can repeat the same line of proof for DSR. Lemma extends immediately 
providing an alternative proof of inequalities (|^. 

Extension of Theorem]^ (i.e., proving that equal vertex weights in VWWI make 
the best estimate of DSR) can be useful in general, but in the context of the present 
paper is a side line of the analysis and can be skipped. 

5 Quality of Lower Bound 

In the previous sections we presented the lower bound TLB{d) and the upper bound 
TUB{d) for the minimum terminal distance spectral radius of a tree with the 
given degree sequence d. We have shown that the terminal distance spectral radius 
TDSR{T**) of such an extremal tree T** G T(d) lies between the average terminal 
distance row sum and the terminal distance spectral radius of the BFS-tiee (which 
can be thought as a “maximally balanced tree” for the degree sequence d) . 

When the average row sum (or, equivalently, the value of ^TWI{BFS(d))) is 
equal to TDSR{BFS{d)), Conjectureholds, and the BFS-tiee minimizes TDSR. 
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But, by the Perron-Frobenius theorem, it is true only when BFS-tiee is terminal 
distance regular (i.e., all row sums of the terminal distance matrix are equal to each 
other), which is not a typical case. 

Nevertheless, below we, to some extent, justify that even when the lower bound 
is not attained, its relative error 


TErr{d) := 


TUB{d) - TLB{d) 
TUB{d) 


( 8 ) 


is small (namely, that it does not exceed 3%). 

Firstly, we provide some computational evidence. In Fig. [^the relative error (|^ 
is presented for all possible degree sequences with the tree order N not exceeding 22. 



Fig. 2. Relative error of the lower bound for TDSR vs the tree order N. 


We see that the relative error never exceeds 3% with most degree sequences having 
error less than 1%. Analysis of “extremal” degree sequences (those constituting the 
upper envelope of the point set in the hgure), shows that all they have the form 


d(a, b) := (1,..., 1, 2,..., 2, a + 6), where a, 6 G N (see Fig. 


a+6 


BFS-tiee for the degree sequence d(3,4)). This observation 
conjecture. 


for the example of the 
eads us to the following 
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Fig. 3. i?FS'-tree for d(3,4). 


Conjecture 3. For every tree degree sequence d there exist such a, 6 G N that 
TErr{BFS{d)) < TErr{BFS{d{a,h))). 

In the following lemma we obtain a closed-form expression for TDSR of the BES- 
tree for the degree sequence d(a, b). 

Lemma 2. 

TDSR{BES{d{a, b))) = 2a + b-3 + + b'^ + 5ab -4a + 2b + l. (9) 

Proof. Denote with R the identity matrix of order n and let Emn stand for the all-ones 
m X n matrix. Consider the terminal distance matrix RD of the tree BES{d{a, b)): 



( 10 ) 


Consider a Perron vector u{RD) = x = (xi, ...,Xa+b)- Since any two of the hrst 
a pendent vertices in BES{d{a,b)) are swapped with no impact on RD, we have 
xi = ... = Xa ='■ OL. Also, the last b vertices in BES{d{a, b)) are similar to each other 
and, hence, Xa+i = ... = Xa+b =■ R 

As, by definition, Ai(i?D)x = RDx, we obtain the following system of equations: 


Xi{RD)a = 4(a — 1)q: -|- 3b(3 
XiIrD)( 3 = 3aa -|- 2(6 — 1)/? 


After elimination of a and jd we find Xi{RD) as the positive root of the square 


equation — 2(2a + b — 3)A — 8(a -|- 6 — 1) — a6 = 0, which immediately gives ([^. □ 


To write down TErr{d{a,b)), evaluate TLB{d{a,b)) as the average row sum of 


RD: 


TLB{d{a,b)) 


4a(a — 1) -|- 2 ■ 3ab -|- 26(6 — 1) 

fl -|- 6 



( 11 ) 


Finally, put the following estimate of TErr{-), which is true if Conjecture]^ holds. 
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Theorem 3. If Conjecture holds, then for any tree degree sequence d 

TErr(d) < -< 0.03. 

3^2 + 4 

Proof. Using Lemma and expression (11) write 


TUB{d{a,b)) 


3 I /i I 0,6 I 26—4a+l 

■r \/ -L + iTTTi: -r 


2a+6 


2 ci -|-6 2 ti -|-6 


Introduce the new variable, n := a + b, and exclude b from (12): 

1-4 

TErr{d{a, n — a)) = 1 — 2- 


a(n-a) , 2n+l-6a 


1 _ o I /1 I _i_ _ 

n-\-a ' \l ' {n+a)2 (n+a)^ 


( 12 ) 


(13) 


Find the maximum by a of the right-hand side of (13) assuming n constant. From 


the hrst-order conditions, the extremal a satishes the equation 


6a/ 3an -|- — 6a -|- 2n -|- 1 -|- — 3an -|- 6a = lOn -|- 2. 


Omitting routine calculations, we conclude that for n > 2 the maximum of (13) is 
attained at a{n) := + 2a/ 2, and 

, . , ^ 8(n-2 + 3^)(l- 4) 

TF'rr(d(a(n), n — a(n))) = 1- - . 

4n - 17 + 6^2 + 3V2n2 - 8n + 6nV2 + 17- 12^2 

Standard analysis shows that TErr{d{a{n),n — a{n))) is a monotone function of 
n. Therefore, for any a, 6 G N 

T£'rr(d(a, 6)) < lim Ti7rr(d(a(n), n — a(n))) = —^-. 

n^+oo 3^2 4 


□ 


6 Conclusion 

Above we suggested a lower bound for the terminal distance spectral radius of a tree 
with the given degree sequence and showed it to be within 3% from the terminal 
distance spectral radius of the BES-tree. This means that Conjecture]^ which says 
that the BES-tree has minimal TDSR among all tree with the given degree sequence, 
is valid at least up to 3% (and typically even more precisely, as shown in Fig. [^. 
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However, our proof of this 3% error relies on the Conjecture which guesses the 
shape of trees, which give the maximum error. Although being pretty natural and 
well-grounded numerically, this conjecture is still an open issue in our analysis. 

It is of interest to discuss how much the presented technique can be applied to 
estimate the precision of Conjecture which states a similar minimum property of 
the BFS-tree with respect to the distance spectral radius {DSR). 

In Fig. I^the relative error 


Err{d) : = 


UB{d{a,b)) - LB{d{a,b)) 


UB{d{a,b)) 

where LB and UB are dehned in (|^, is presented for all trees of order N < 23. We 
see that the error of LB is, in average, twice as big as that of TLB (see Fig. [^. From 
Fig. [^one might conjecture that Err does not exceed 6%, but careful analysis, similar 
to that performed in Lemma needs to be performed to justify this conjecture. 

An observation, which may help, is that degree sequences giving the maximum 
Err (and, hence, forming the upper envelope of the point set in Fig. have the form 
(1,..., 1,2,..., 2, d). They correspond to, the so-called, starlike trees, whose distance 


d times 


spectral properties are studied in detail by Stevanovic and Ilic 17 



Fig. 4. Relative error Err of the lower bound for DSR vs the tree order N. 
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7 Appendix. Proof of Theorem [2 


First we recall the notions of a directed tree and snbordinate gronp weight from [^. 

Definition 5. A vertex-weighted directed tree is a connected directed graph where 
each vertex except the root has the sole outbound arc, the root has no outbound arcs, 
and each vertex has non-negative weight assigned to it. 


An arbitrary vertex-weighted tree T consisting of more than two vertices can be 
transformed into a directed tree by choosing an internal vertex r G M{T) as a 
root, and replacing all its edges with arcs directed towards the root. Let yV7Z{p., d) 
stand for all directed trees obtained from >VT(/i, d) in this way. 

If at Step i of the Hnffman algorithm (see Section directed arcs towards the 
vertex m* are added instead of undirected edges, the algorithm builds a directed 
Huffman tree. Let TZ'H{pi,d) be the set of directed Huffman trees for a generating 
tuple (/r, d). 


Definition 6. For an arbitrary vertex v G V (T) of the directed tree T define its 
subordinate group griv) F V{T) as the set of vertices having the directed path to 
the vertex v in the tree T {the vertex v itself belongs to griv))- The weight friv) 
of the subordinate group griv) is defined as the total vertex weight of the group: 
Mv) ■■= J2uGgTiv)Mu)- 

Definition 7. For a directed tree T G W7^(/i, d) with root r define a vector f(T) : = 
(/r("i))mgA^(r)\{r} of subordinate groups’ weights. 

If some tree T G WT{pi, d) is transformed into a directed tree Tr G WlZ{pi, d) by 
choosing a root r, its Wiener index can be written as [I^ : 

VWWI{T) = VWWI{Tr)= fTr{v){fl-fTM)= E xifniv)), 

■yey(r)\{r} •ueV(T)\{r} 

(14) 

where x(a;) := x{p, — x). 

Lemma 3. |[^| If weights are degree-monotone in {fi,d), then for any H G 7l'H{fi,d) 
[vm, v'm' G E{H),m m', fu{v) < fniv')] =» fnim) < /^(m'). (15) 


Definition 8. |15|2^ A non-negative seguence x = {xi,...,Xp), p G N, weakly 


majorizes a non-negative seguence y = {yi, ...,yp) {denote it y ^ or ^ y) if 

k k 


for all k = l,...,p. 

i=l i=l 

//x^ 7 ^ y^, then x is said to strictly weakly majorize y {denote it y -<‘^ ^ or^ >~^ y). 
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We will need the following properties of weak majorization. 


Lemma 4. 1 


Consider a positive number b > 0 and two non-negative sequences, 
X = (xi, ...,Xk,yi, andy = (xi+6, ...,Xk+h,yi-h, ...,yi-b), such thatO <k <1. 

If Xi > yi for f = 1, k, then x y. 


Lemma 5. |15,2^ //x y and x' y', then (x,x') (y,yO; where (x,x') is 
concatenation of sequences x and x'. Moreover, ifx.' y', then (x, x') (y,yO- 


Lemma 6. |15[ , |24| // x{x) is an increasing concave function and {xi,...,Xp) :<^ 
(?/i,I/p), then YfJi=iX{.Xi) > equality is possible only when 

{xi,...,Xpy = (|/i,...,|/p)^. 

Theorem 4. |[^| If weights are degree-monotone in a generating tuple (/r, d) and 
H G IZ'H{yi, d), then for any directed tree T G W7^(/i, d) i{H) f(T). 


Lemma 7. Consider a degree sequence d with n elements being equal to unity and 
a positive sequence pi = (pi,...,/!„) such that p,i > pij for some i and j. Define 
a generating tuple {y',d), which differs from {p,,d) only with i-th and j-th weight 
components, namely, pi[ — e = pit + e = \{pii + pij), where < e < 5 ■.= \{pii — pij) . If 
H G Il'H{p,,d) and H' G lZ'H{p!,d) are some directed Huffman trees, then f(-ff) 

Proof. Suppose vertices u,v E W (/r, d) have weights pj and pij respectively. Let 
(m, Ui,Uk, m) and (n, Vi, vi,m) be the disjoint paths in the directed Huffman tree 
H' G ^^{pifd) from vertices u and v to some vertex m G M{H'), where k,l > 0. 
Since fufu) = pi[,fH'{v) = pit, and pi) > pit, it follows immediately from Lemma 
that k < I and fn'iui) > fwipi) for i = 1, ...k. 

Consider a directed weighted tree T obtained from H' by changing the weights of 
vertices u, v to pii and pij respectively. In the tree T weights of groups subordinated 
to the vertices Ui,...,Uk increase by 6 := 5 — £ > 0 (i.e., friwi) = fwiui) + b,i = 
1,..., k), weights of the groups subordinated to the vertices Vi, ...,vi decrease by b (i.e., 
frivi) = fH'ixi) — b,i = l,...,l), weights of all other vertices (including m) do not 
change. 

Therefore, by Lemma 


y := (/t(mi), •••, friuk), /^(^^l), •••, frivi)) = 


— ifufui) + b ,..., fnfuk) + b, fH'{vi) — b, ...,fH'{vi) — b) 
ifn'iui), ...,fH'{uk), fH'ivi), ...,fHi{vi)) =: x. 


16 







If one denotes with z the sequence of (unchanged) weights of groups subordinated 
to all other internal vertices of T distinct from the root, then, by Lemma § f(T) = 
(y,z) (x,z) = 

It is clear that T G W7^(p, d). Since H G from Theorem]^ we know 

that f(iL) f(T) and, consequently, i{H) □ 

Corollary 1. Under the conditions of Lemma^TVWW I{H) < TVWWI{H'). 

Proof. Let r and r' be roots of trees H and H' respectively. According to 0. 

n 

TVWWI{H) = E xiUk) + Xifniv)), where x{x) = a;(p - x), 

k=l vGM{fi,d)\{r} 

and 

n 

TVWWI{H') = E x(hfc)+ xUH'iv)). 

k=l v&M{ii,d)\{r'} 

From Lemma 0 we know that f(id) Therefore, by Lemma 

Xifniv)) < Y Xifn'iv)). 

v£M{fL,d)\{r} v£M(fi'^d)\{r'} 

Since x(a;) is concave. 


xiUi) + xiUj) < xiUi - d + £) + xiUj + S-e). 


Other vertex weights do not change, so TVWWI{H) < TVWWI{H'). □ 

Lemma 8. Consider a degree sequence d with n elements being equal to unity and a 
vector fi E Sf. If Hi > fij for some i and j, then such a vector u E Sf exists that 
TVWWI{H") < TVWWI{H), where H" E H{yi,d) and H E H{v,d) are Huffman 
trees. 


Proof. Recall a weight vector /i' from Lemma and dehne the vector u E Sf 


as 


= 



+ (p-i + Pj)£, 


k'i = 



(pi + Pj)£, 


k'k — Pki k ij j. 


Interrelation between p, /i', and u is shown in Fig. 

Let finiu) = Pi, ^h{v) = Hj for some u,v E V{H) and consider a vertex-weighted 
tree T obtained from H by changing weights of u and v to p' and p)- respectively. 
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From Fig. 2 is is clear that z/j > /i', Vj > jj,'-. Since RD{-) is off-diagonal positive, 
TVWWI{-) is strictly monotone in weights, so TVWWI{T) < TVWWI{H). Con¬ 
sider a tree H' G By Note || we have TVWWI{H') < TVWWI{T). From 

Corollary we also know that TVWWI{H”) < TVWW Summarizing these 
inequalities obtain the statement of the lemma. □ 



Fig. 5. Interrelation between /r, /i',and z/. 

Lemma 1^ says that when positive weights for pendent vertices are picked from the 
unit sphere, TVWWI{-) of a tree with the degree sequence achieves its maximum 
when all vertex weights are equal to each other, i.e., pj = I/a/u for all i = 1, 

The vector e := of equal weights is the only positive vector on the 

unit sphere, to which Lemma cannot be applied, and, since for an arbitrary graph 
G the identity holds e^RD{G)e = l"'"RD{G)l/n = 2TWI{G)/n, the the statement 
of Theorem follows immediately. 
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